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ABSTRACT 

We present a detailed study of the small frequency separations as diagnostics of the 
mass of the convective core and evolutionary stage of solar-type stars. We demonstrate 
how the small separations can be combined to provide sensitive tests for the presence 
of convective overshoot at the edge of the core. These studies are focused on low 
degree oscillation modes, the only modes expected to be detected in distant stars. 
Using simulated data with realistic errors, we find that the mass of the convective 
core can be estimated to within 5% if the total stellar mass is known. Systematic 
errors arising due to uncertainty in the mass could be up to 20%. The evolutionary 
stage of the star, determined in terms of the central hydrogen abundance using our 
proposed technique, however, is much less sensitive to the mass estimate. 
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1 INTRODUCTION 



The theory of stellar structure predicts that the central re- 
gion of massive stars (M ^ I.IMq) are convective, rather 
than radiative as is the case for the Sun. Since convective 
flows imply chemical mixing, the evolution of these stars is 
severely influenced by the presence, and the extent, of the 
convective core. The process of convective transport also in- 
volves the question of overshoot which implies the extension 
of the convective region beyond the classical border of con- 
vective stability. This causes more hydrogen to be supplied 
to the nuclear reactions going on in the core, thus effectively 
increasing the lifetime of a star on the main sequence. The 
present theories of convection are inadequate to provide a 
definite measure of the extent of overshoot and it is one of 
the free parameters in a stellar model. The uncertainty in 
the mass of the convective core, or its extension due to over- 
shoot, directly translates into an error in the determination 
of the age of the star. This in turn would affect other stud- 
ies based on stellar evolution like population synthesis in 
galaxies etc. 

The size of the convective core does not remain con- 
stant, and changes as the star evolves. The variation in the 
core size of a I.SMq star is illustrated in Fig. ^through the 
change in the hydrogen profile as a function of fractional 
mass. For stars of mass 1.1 M/Mq ^ 2.0 the core in- 
creases in size during the initial stages of evolution before 
beginning to shrink later. Furthermore, stellar models indi- 
cate that the extent and evolution of the convective core is 
quite different depending on whether or not convective over- 
shoot is present (Fig. It is, therefore, clear that a good 
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test of the presence and extent of convective layers in the 
central region of massive stars will go a long way towards 
our understanding of stellar structure and evolution. While 
there is no direct method to look inside a star, the recent 
advances in asteroseismology does offer an opportunity to 
probe the innermost layers of stars. The frequencies of os- 
cillation of a star depend on its global properties like the 
mass, radius etc, as well as the detailed internal structure, 
especially the location of layers of rapid change in density or 
chemical composition, and hence can be used to probe the 
inner layers. 

Space-based a steroseismic missions suc h as MOST 
(Walker et al.ll2003l) and CoRoT teaglinll2003l) are expected 
to measure the oscillation frequencies of several stars with 
sufficient accuracy to enable a more detailed study of stel- 
lar structure than has been possible so far. Ground-based 
observations have already started measuring the frequencies 
of many stars despite the difficulties of such measurements. 
Unfortunately, none of these observations can or will be able 
to observe the intermediate degree modes of oscillation that 
have been so useful in solar studies. These missions can at 
most, hope to determine the frequencies of modes with de- 
gree £ Si 3. However, we are fortunate in the fact that the 
modes of low degrees (I = 0-3) penetrate the central regions 
of a star and hence carry information about its deepest lay- 
ers. We focus our study only on these modes. 

It has been suggested that the oscillatory signal in 
the frequencies due to sharp features such as the bound- 
aries of convective regions or ionisation zones can be 
used to estimate the location of these layers inside a 
star iMonteiro. Christensen-D alsgaar d fc Thompson! Il998l . 
2000). While this technique appears to succeed for the base 
of the outer convective envelope or the second helium ionisa- 
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Figure 1. The hydrogen profile near the boundary of the convec- 
tive core is shown for a 1.3Mq model with overshoot at different 
stages of evolution. The fiat portion of the curve indicates the 
core which is fully mixed due to convection. It is clear that the 
size of the core increases during the first part of evolution before 
shrinking in later stages. Diffusion of hydrogen leaves a smoothed 
out step-like signature in the hydrogen profile as the core recedes 
with increasing age of the star. 



tion z one jBallot, Turck-Chieze. fc Garciall2004l ; iBasu et alJ 
l2004ft . lMazumdar fc Antial feOOll) have shown that it would 
not work for the convective core due to a problem similar 
to aliasing in Fourier analysis. In this paper we show that 
combinations of frequency separations can provide a better 
knowledge of the core. 

It is well known that the large scale properties like the 
mass and radius act mostly to scale the frequencies of nearly 
all modes of oscillation. Such scaling effects are apparent 
in the large separation, i.e., the difference in frequency of 
modes of successive order but of a particular degree. The 
large separations are also quite sensitive to the structure of 
the outer layers of a star, which are among the most uncer- 
tain aspects of theoretical stellar models. On the other hand, 



0.10 



0.08 



0.06 



0.04 



0.02 



0.00 




0.6 



0.4 



0.2 



Figure 2. The evolution of the fractional mass of the convective 
core, m COTe /M is shown for models of mass 1.2, 1.3 and 1.4Mg . 
The black lines represent models without convective core over- 
shoot and the other curves show models in which core overshoot 
(0.2Hp) has been incorporated. 



the so-called small separations, i.e., the small differences of 
frequency between modes of nearly same order but different 
degrees, are more affected by the deeper layers of a star. 
Therefore, the small separations are useful in investigating 
the central features, especially the extent of convective cores 
in massive stars. 

Indeed, many efforts have been made to link the 
small frequency separations to the convective core 
in massive stars (for theoretical studies see, e.g., 
[Dziembowski fc Pamvatnvkh| | l99l| j [Roxburgh fc Vorontsovl 
hoOll) . Using a Cen models, buenther fc Demarcate! ll200m 
had suggested to use the small separations to constrain the 
age of the st ars. After oscillations were indeed detected 
in a Cen AB (|Biiuchv_J^XJamer||2QQ3; ICarri er fc BourbarJ 
1200.4 iBedding et alJ 12004 iKieldsen et al.ll2005T l. different 
authors have made use of the small separations to con- 
strain the parameter s , especially the age, for these stars 
jThevenin et alJl2002l iThoul et al.ll2005l:jEggenberger et"aH 
2004'!. Recentlv. lRoxburgh fc VorontsovT(l2003t) have demon- 
strated that the ratio of the small and large separations 
serve as a better diagnostic of the i nnermost layers of a 
star t han the small separations alone. iMiglio fc Montalbanl 
(2005) have used this technique to refine the seismic 
models for a Cen AB. Similarly, the small separations 
have proved to be usefu l for the seismic modelling of 
Procy on A as well ijEp jpjenberpje^ Carr ier, fc B ouchvl 
[ 20051: iProvost. Berthomieu. fc Martid ~ l200rj) . 

IStraka. Demarau^^htentff^'il I2OO5I) have provided 
upper limits to the convective overshoot in Procyon A 
models using the frequency separations. In these studies, 
the frequency separations have been mostly used to find the 
best possible model that fits the observed data for a given 
star. In this paper we aim to present a general technique 
for probing specific properties of the central regions of a 
star, like the mass of the convective core, and the state of 
evolution of the star, using suitable combinations of the 
frequency separations. 
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While the large and small separation corresponding to 
each radial order encodes information about a different layer 
inside the star, and is, therefore, an independent diagnostic 
of the structure, there are also advantages in considering 
suitable average values of these separations. Firstly, the av- 
eraging over several radial orders removes the small-scale 
variations due to subtle differences in the internal struc- 
ture, and enhances the effects of the more significant fea- 
tures of the stellar interior. Second: it is often not possible 
from observations to determine the frequency of modes of 
all radial orders for the same degree, and an average is of- 
ten the best indication of the trend of frequencies. Third: 
while individual frequencies, and hence the separations are 
often susceptible to large observational errors, the average 
values of the lar ge and small sepa rations turn out to be ro- 
bust quantities (lMazumdaill2005T) . The diagnostic power of 
the average separations can be util ised very well through 
the C hristensen-Dalsgaard diagram dChristensen-Dalsgaardl 
1988) to determine the mass and age of a star. In this work, 
we show that they can also be used for more detailed study 
of the stellar interior and present a practical method involv- 
ing the average values of small and large separations and 
combinations thereof to estimate the age of a star and the 
mass of the convective core. 

In the next section we define new diagnostic quanti- 
ties to probe the stellar core and outline our technique and 
describe the stellar models used in this work. Our results 
along with error estimates are presented in Section [3] with 
the concluding discussion in Section 2] 



2 THE TECHNIQUE 

The average large separation of radial modes is defined as 

A = Hn + l,O)-Kn,0», (1) 

where the angular brackets imply averaging over a range of 
radial order, n, appropriately chosen. 

The small separation between modes of degree £ and 
£ + 2 is defined as 



See +2 (n) = v{n, I) - v{n - 1,1 + 2) 



(2) 



While comparing small separations between different pairs 
of modes (e.g., (£ — 0,2) and {£ — 1,3)) it is convenient 
to scale the differences to eliminate the factor of £{£ + 1) in 
the frequencies. Thus one may define the scaled small sepa- 
rations as (cf. [Chns^enseii-Da lsgaard fc Berthomieulll99lt 
iRoxburgh fc Vororitsovll200^1) 



n l„\ S "+2( n ) v{n,£) - v{ 
+ K 1 41 + 6 A" 



n- l,£ + 2) 



(•3) 



In the rest of the paper, we will refer to this scaled defi- 
nition, Dee+2, simply as the "small separation". The small 
separation between modes of degree £ = and 2 are des- 
igna ted as Dq 2 , and those correspo nding to £ = 1 and 3 as 
-Dm. lRoxburgh fc Vorontsovl (l200ll) have demonstrated that 
these two sets of small separations encode similar, but some- 
what different information about the central layers. This is 
explained by the fact that modes of different degree have 
different inner turning points, and therefore probe different 
layers. Thus, it might be useful to compare the two small 
separations, Do 2 (n) and Di3(n), of the same radial order. 



(4) 
(5) 



In order to exploit the fine differences between Do2(n) 
and Dis(n), for each model, we define the following combi- 
nations of the small separations: 

e( ni ,n 2 ) = ^ {D 13 (n) - D 02 (n)) n=niton2 

Vim,**) = (^) 2 <(A,(») - D 02 (n)f)_ niton2 

where the arguments n\ and n 2 indicate the limits of radial 
orders between which the averaging is carried out. Aq is the 
average large separation of the Sun (A0 ~ 135/iHz). 

The properties of the quantities 8 and rj, and their sen- 
sitivity to different stellar parameters, would partly depend 
on the range of the radial order chosen for averaging. We 
deliberately choose relatively higher order modes (n ^ 10) 
for averaging to avoid contamination from mixed modes. 
In practice, however, the choice of the range for averaging 
would be primarily governed by the availability of data in 
different frequency ranges of the observed spectrum. For a 
complete set of frequencies of a solar-type star, we find that 
the range of n = 15 to 28 for 6 and n = 10 to 28 for rj make 
these quantities ideally sensitive to stellar properties such 
as the convective core mass, m core or the central hydrogen 
abundance, X c . In the rest of the paper, we shall assume 
these to be the respective ranges of averaging for and rj, 
unless the range is stated explicitly. 

We include the scaling factors Aq /A in the definitions 
of 8 and rj to reduce homology effects while comparing stars 
of different mass and radius. This factor essentially filters 
out the global scaling of the frequencies due to differences in 
the mean density of different stars, and enhances the effects 
of the structural details which we want to study. 

We study the variation of these new diagnostic quanti- 
ties, 8 and rj with the age and evolution of the convective 
core in stars using a large number of stellar models. We re- 
strict ourselves to stars with masses just a little higher than 
the mass of the Sun, namely, the range 1.2-1.4Mq . We study 
stellar models evolved to different ages on the main sequence 
for masses 1.2, 1.3 and IAMq. Ideally, the boundary of the 
growing convective core should be tracked by adjusting the 
mesh at each evolutionary step. However, since the current 
evolutionary code does not have that capability, the edge of 
the core is determined through the Schwarzschild criterion 
being tested over a very fine mesh in the central regions of 
the star. To compare with stars without convection in their 
cores, we also constructed a sequence of IMq models. We 
also vary the convective overshoot in the core, which is mea- 
sured in units of the local pressure scale height, Hp. We use 
two sets of models with the overshoot parameter d ov = 
and d ov = 0.2, where the amount of overshoot at the edge of 
the convective core is given by d ov Hp. For very small con- 
vective cores, the overshoot is limited to be a fraction of the 
core size instead of Hp. However, for most of the models 
that we consider, the size of the core is typically at least an 
order of magnitude larger than the pressure scale height at 
the boundary, so that we can simply use d ov Hp as a measure 
of the overshoot. 

The models used in this work were constructed us- 
ing YREC, the Yale Rotating Evolution C ode in its non- 
rotating configuration jGuenther et alJll992ft. These models 
use the OPAL eq uation of state (iRqgers fc Navfonovll2002ri . 
OPAL opacities lllglesias fc Rogers!E^9^) . low temperature 
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Figure 3. The evolution of the average small separation {D02) 
(upper panel) and the new diagnostic r\ (lower panel) as a func- 
tion of the effective temperature are shown. The black lines rep- 
resent models without overshoot while the other lines represent 
models with overshoot (d ov = 0.2). For each mass the curves 
for models with and without overshoot diverge from the zero age 
main sequence. The errorbars represent la errors in the respective 
quantities for relative frequency errors of 10 -4 , obtained through 
simulations described in Section kj.4,11 The significant kink in the 
curves, especially for overshoot models, correspond to the age of 
the star when the convective core stops expanding and begins 
shrinking. 



opacities of lAlexander fc Ferguson! lll994h and nuclear reac- 
tion rates as used bv TBahTall^^Pinsonneaultl (^992). The 
models take into account diffusio n of helium and heavy 
eleme nts, using the prescription of lThoul. Bahcall fc Loea 
( 1994). The initial chemical composition of the models was 
fixed at (Y = 0.26, Z = 0.022). 



3 RESULTS 

3.1 Testing for core overshoot 

As is clear from Fig. |5] the evolution of the central con- 
vective region of a star with age is significantly affected by 
the presence of overshoot. Indeed, we find that this impor- 
tant difference is reflected in the oscillation frequencies and 
consequently, the small separations. The presence of convec- 
tive overshoot, therefore, might be tested through the small 



separations themselves. iDi Mauro et alJ i200ot) have made a 
detailed investigation on how the frequency separations can 
be used to distinguish between models with and withou t 
overshoot for a post- main sequence star. iMazumdarl i2005t) . 
however, has shown that during the main sequence phase 
the direct effect of overshoot on the small separations D02 
becomes noticeable only in the later stages of evolution. We 
investigate whether the average small separations, (-D02), or 
the quantity r\ can provide any clue to the presence of over- 
shoot in a given star. 

In Fig. |HJ we show how the average small separations 
{D02) and the quantity r\ vary as functions of the effective 
temperature as the star evolves on the main sequence, for 
models with and without core overshoot. The errorbars rep- 
resent the la errors in the ordinates for relative frequency 
errors of 10 -4 . These errorbars were obtained from Monte 
Carlo simulations with theoretical frequencies, described in 
detail in Section 13.4.11 It is evident that even if the mass 
of a star is known, (D02) is barely sufficient to distinguish 
between models with and without overshoot at a given T e g. 
The tracks for rj, on the other hand, are well separated, 
with the overshoot models producing much higher values of 
rj at a given T e g , irrespective of mass. While the {-D02) val- 
ues for models with and without overshoot are nearly the 
same, especially near the zero age main sequence (ZAMS), 
the rj values of overshoot models diverge very rapidly be- 
yond the ZAMS from those of non-overshoot models. Thus, 
if the location of a star on the HR diagram is known, the 
quantity r\ can be used to test the presence of overshoot in 
the convective core. For large errors in frequency the errors 
in 77 increase almost linearly (see Table 0; even for relative 
errors 5 x 10 -4 the separation between the tracks is large 
enough to be distinctive to prove useful. 



3.2 Estimating the convective core mass, m coro 

As evident from Figs.0and[5] the fractional mass of the con- 
vective core of a star of total mass of about 1.2-1.4M0 does 
not remain constant throughout its main sequence evolution. 
In particular, the convective core grows during the early part 
of evolution and then begins to shrink at later stages. The 
exact evolution of the core depends somewhat on the total 
mass and quite strongly on whether or not convective over- 
shoot occurs. In either case, the relationship between the 
fractional mass of the core and the age of the star has op- 
posite signs in two parts of evolution on the main sequence. 
Since the mass of the convective core is not monotonic with 
the age, it is difficult to associate the core mass directly with 
the small separations (and hence with 9 or 77), which vary 
generally monotonically with age. However, we find that if 
we scale the mass of the core with the factor (Aq /A) 2 which 
is related to the mean density of the star, the scaled mass is 
nearly monotonic with the age of the star. Of course, even 
this scaled mass of the core eventually decreases with in- 
creasing age towards the end of the main sequence when the 
core shrinks to a very small size. While theoretically this 
might appear to be somewhat ad hoc, it provides a practical 
way to overcome the difficult dichotomy of the core mass. 
Observationally, it is a perfectly suitable method since the 
average large separation can usually be measured to a fair 
degree of accuracy. If the scaled core mass can be determined 



5 



0.2 - 



0.1 



8 0.0 





< 0.4 




0.2 



2.0 


1.5 


1.0 


0.5 


i i i 1 


i i i 1 

1.2 M 

1.3 Mg, 

1.4 M e 


1 1 1 | 1 1 

y 

/* 

/ 


i i 1 i i 

/ 






/ / 


















= o.2 : 


i 1 


i i i 1 


i i i i i i 


i i i i i 


2.0 


1.5 
(D 02 


1.0 
) (MHz) 


0.5 



0.2 - 



0.1 



8 o.o 
£ 

<r 

\ 



< 0.4 




d„, = 0.0 



J i i i I i i i I i i i L 



J I I i L 



0.02 0.04 0.06 0.08 0.10 



0.1 - 



T 



1 — 1 — 1 — 1 — I — 
1.2 M g 

1.3 M G 

1.4 M B 




d... = 0.2 



J I I I I I I I I I I I L 



0.2 0.4 0.6 

77(10,28) (mHz 2 ) 



Figure 4. The scaled mass of the convective core, 
(A0/A) 2 m C orc/Af is plotted as a function of the average small 
separation, {D02), for models without overshoot (upper panel) 
and with overshoot (lower panel). The dotted line shows the (D02) 
values for a IMq model, for which the core is not convective. The 
errorbars represent lc errors in (D02) for relative frequency errors 
of 10 — 4 , obtained through simulations described in Section l3.4.1l 
The significant kink in the curves, especially for overshoot mod- 
els, correspond to the age of the star when the convective core 
stops expanding and begins shrinking. 



Figure 5. The scaled mass of the convective core, 
(A0/A) 2 m C orc/Af is plotted as a function of the diagnostic, rj, 
for models without overshoot (upper panel) and with overshoot 
(lower panel). The dotted line shows the n values for a IMq 
model, for which the core is not convective. The errorbars repre- 
sent 1(7 errors in n for relative frequency errors of 10 -4 , obtained 
through simulations described in Section 13.4.11 The significant 
kink in the curves, especially for overshoot models, correspond to 
the age of the star when the convective core stops expanding and 
begins shrinking. 



using a seismic diagnostic, the actual mass of the core can 
easily be extracted through the average large separation. 

We have also found from Fig. [3] that the value of r\ is 
strongly dependent on the presence of overshoot. This im- 
plies that we need separate calibration curves for models 
with and without core overshoot. Since Fig. ^itself provides 
a means of testing for the presence of overshoot, from a prac- 
tical viewpoint, this does not pose a problem while dealing 
with a given star. We use multiple models at different ages 
of 1.2, 1.3 and IAMq both with and without overshoot to 
produce the needed calibration curves. 

We first investigate the variation of the average small 
separation, (D02) with the scaled fractional core mass, 
(A /A) 2 m coro /M. In Fig.[| we show the scaled core mass as 
function of (-D02), for stellar models of different mass, evolv- 
ing through the main sequence. Every curve in this figure is 
an evolutionary track, which serves as a calibration curve to 
determine (A0/A) 2 m CO re/A^ for a given target. The figure 
(and similarly Figs [5] |S| and |7J have been drawn such that 



the sense of evolution is from the left to the right along the 
x-axis. 

This figure illustrates the fact that although models 
with the same core mass have largely different values of 
the small separation, the scaled core mass is almost mono- 
tonic with {Dq2}- It is worth mentioning that the kink in the 
curves shown in the lower panel of Fig. [I] is not an artifact 
of the stellar models. It is, in fact, intrinsic to the evolution 
of the star, and corresponds to the age when the increase of 
the core size stops and the core begins to shrink (cf. Fig.|2J. 
Such a kink is also seen in the Christensen-Dalsgaard dia- 
grams jMazumdarll2005D . The presence of the kink implies 
that a careful calibration needs to be carried out - perhaps 
two different calibrations for two segments of the evolution- 
ary track. For comparison, we have also plotted the (D02) 
values for a IMq model, which does not have a convective 
core at all. Every horizontal errorbar indicates the la er- 
ror in {D02) due to relative frequency errors of 10" 4 , which 
would translate into a random error in the determination of 



6 A. Mazumdar et al. 




Figure 6. The central hydrogen abundance, X c is plotted as a 
function of the diagnostic, rj, for models without overshoot (up- 
per panel) and with overshoot (lower panel). The errorbars repre- 
sent la errors in n for relative frequency errors of 10 — 4 , obtained 
through simulations described in Section 13,4.11 The significant 
kink in the curves, especially for overshoot models, correspond to 
the age of the star when the convective core stops expanding and 
begins shrinking. For the no-overshoot models rj is poorly sensi- 
tive to X c at early stages of evolution and hence it will not be 
possible to determine X c for low values of r). 

(A0/A) 2 m CO rc/M through the local slope of the calibration 
curves. Given the wide separation of the curves correspond- 
ing to different masses, it is clear that such random error in 
estimating the core mass due to errors in frequencies is much 
smaller than the uncertainty due to the stellar mass. There- 
fore, to use this technique, one needs to have a fair idea of 
the mass of the star. We shall investigate both the random 
and systematic errors involved in detail in Sections 13.4.11 
and 13.4.21 respectively. 

We can also use the quantity r\ as an additional di- 
agnostic for the convective core. The variation of 77 with 
(A /A) 2 m corc /M is illustrated in Fig. El We find that r) is 
more sensitive to the core mass in the presence of overshoot. 
Actually, the use of 77 as a diagnostic is restricted only to 
more evolved models in the absence of overshoot, especially 
for the 1.4Mq models. The presence of the kink again im- 
plies that for stars near that particular evolutionary stage 
the error in the core mass estimate will be larger. We also 
note that unlike {D02), the range of values of 77 for a star 



Figure 7. The central hydrogen abundance, X c is plotted as a 
function of the diagnostic, 9, for models without overshoot (up- 
per panel) and with overshoot (lower panel). The errorbars rep- 
resent la errors in 9 for relative frequency errors of 10" 4 , ob- 
tained through simulations described in Section 13.4.11 It turns 
out that the correlation of 9 with X c becomes poor beyond a cer- 
tain evolutionary stage (depending on the mass) and the curves 
are truncated only to show the part that would be usable for the 
determination of X c . 

without a convective core, the IMq model, is restricted to 
much smaller values than that for stars with a convective 
core. This is one advantage that 77 has, as a diagnostic of 
m corc , over (D02) (cf. Fig. 0J. 

3.3 Estimating X c 

The central parts of a star change more rapidly than the 
outer layers with age. Thus, it is not surprising that the 
small separations, or co mbinations thereof, can be us ed as 
indicators of stellar age ( Christcnsen-DalsEaard 1988). We 
have investigated how sensitive the quantities 77 and 9 are to 
the stellar age, as measured in terms of the central hydrogen 
abundance, X c . Figs.|S|and|7|show the calibration curves. As 
with the core mass, we find that in the absence of overshoot 
77 is not very sensitive to X c either. For overshoot models, 
however, 77 may be used as an indicator of X c throughout the 
main sequence phase, provided the mass is independently 
known. The quantity, 6, on the other hand, serves as a good 
diagnostic for X c in the early stages of evolution. Beyond a 
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Figure 8. The mass of the convective core, as determined 
through the calibration of (A0/A) 2 m CO ro/-M vs (D02) (Fig.|3J, is 
plotted against the model value. Results are shown for test models 
both without (upper panel) and with (lower panel) overshoot. In 
each case, the calibration curve corresponding to the same mass 
and overshoot was used. The dotted lines show the values which 
would be obtained for a perfect match with the model. The er- 
rorbars show the 3c error estimates for relative frequency errors 
of 10" 4 . 



Figure 9. The mass of the convective core, as determined 
through the calibration of (Aq / A) 2 m core /M vs n (Fig. [SJ, is 
plotted against the model value. Results are shown for test mod- 
els both without (upper panel) and with (lower panel) overshoot. 
In each case, the calibration curve corresponding to the same mass 
and overshoot was used. The dotted lines show the values which 
would be obtained for a perfect match with the model. The er- 
rorbars show the 3c error estimates for relative frequency errors 
of 10" 4 . 



certain age, which depends slightly on the mass, 6 no longer 
follows a regular pattern with X c . At younger ages, however, 
9 is less sensitive to the stellar mass, which would make it 
an excellent indicator of the evolutionary stage if the mass 
is not well-constrained. The quantities and r\ also serve 
as complementary diagnostics of X c at different stages of 
evolution. 

3.4 Errors in estimation of m CO r C and X c 
3.4-1 Random errors 

Having constructed calibration curves for the core mass and 
stellar age, the next step is to check whether we can use 
them for real data with errors. Since at present we do not 
have enough observational data to test these methods, we 
undertook a Monte-Carlo simulation where random errors 
with a Gaussian distribution that has a specified standard 
deviation were added to the exact frequencies of a given test 
model. All the necessary separations were then calculated 



with these error-added frequencies and the desired quantity 
(wicore or A c ) was estimated using the calibration curves. 
Finally, the averages of the estimates along with errorbars 
were obtained from 100 such realisations. Typically we used 
five test models at different evolutionary stages for each mass 
to check the viability of the method. For each test model, it 
is assumed that the mass of the target star is known and that 
the question of overshoot has been settled independently - 
through a test like Fig. |H] or by some other constraint - 
so that we can use a calibration curve constructed for the 
same mass and overshoot. The systematic errors arising due 
to uncertainties in these respects will be treated separately 
in the next section. 

Figures 151 to 1 1 1 1 show the results of these tests. In each 
figure we compare the original input model value of m CO r C 
and X c against their estimated values. The errorbars indi- 
cate 3cr errors obtained assuming a constant relative error 
of 1 part in 10 4 in frequencies. Such errors in frequencies 
are in line with the expectations from the C0R0T mission. 
Even recent ground-based observations have reached almost 
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Table 1. Results of Monte Carlo simulations of determination of m COTC and X c for 5 test models at different ages with mass M = 1.3Mq 
and overshoot d ov = 0.2 are shown. For each model the top row shows the original input model values. The three subsequent rows show 
the results for different quantities when relative errors in frequency of 1 X 10 — 4 , 5 X 10~ 4 and 1 X 10 -3 respectively are added to the 
model frequencies. In each case, the results represent the average and 1<t errors for 100 independent realisations of the simulated data. 
The blank entries indicate that the technique fails to extract the desired quantity either due to lack of sensitivity of the diagnostic, or 
too high an error in frequencies to use the calibration curves. 



Age 
(Gyr) 


A 

(MHz) 


(D02> 
(MHz) 


e 

(MHz) 


n 

(/■Hz 2 ) 


m coro /M 
using using 
(D a) V 


x c 

using using 
V 8 


1.300 


108.862 


1.689 


0.093 


0.020 


0.083 


0.604 




108.863±0.017 
108.869±0.086 
108.876±0.173 


1.689±0.014 
1.688±0.067 
1.687±0.134 


0.094±0.025 
0.099±0.126 
0.106±0.252 


0.026±0.005 
().171±0.052 
0.620±0.195 


0.084±0.001 0.088±0.003 
0.084±0.003 0.143±0.018 
0.084±0.007 


0.585±0.012 0.602±0.014 
0.394±0.063 0.599±0.070 
0.578±0.096 


1.885 


102.054 


1.477 


0.212 


0.042 


0.086 


0.544 




102.055±0.016 
102.060±0.081 
102.067±0.162 


1.477±0.012 
1.475±0.062 
1.474±0.124 


0.213±0.025 
0.216±0.126 
0.221±0.251 


0.048±0.009 
0.195±0.070 
0.650±0.234 


0.084±0.001 0.087±0.004 
0.084±0.004 0.131±0.017 
0.085±0.007 


0.540±0.018 0.536±0.014 
0.370±0.074 0.539±0.066 
0.576±0.092 


3.000 


88.636 


0.957 


0.408 


0.148 


0.092 


0.411 




88.637±0.014 
88.642±0.070 
88.648±0.141 


0.957±0.011 
0.956±0.055 
0.955±0.109 


0.409±0.025 
0.414±0.126 
0.421±0.252 


0.156±0.017 
0.304±0.098 
0.759±0.258 


0.088±0.001 0.090±0.008 
0.090±0.005 0.112±0.011 
0.091±0.009 


0.420±0.037 0.435±0.010 
0.281±0.080 


3.855 


78.129 


0.636 


0.505 


0.275 


0.089 


0.280 




78.130±0.013 
78.134±0.062 
78.139±0.124 


0.636±0.010 
0.635±0.048 
0.635±0.097 


0.506±0.026 
0.510±0.128 
0.515±0.255 


0.282±0.025 
0.433±0.133 
0.897±0.322 


0.088±0.001 0.088±0.002 
0.088±0.002 0.092±0.006 
0.089±0.004 


0.286±0.022 
0.223±0.068 


4.400 


71.480 


0.428 


0.550 


0.407 


0.080 


0.173 




71.481±0.011 
71.485±0.057 
71.490±0.114 


0.428±0.009 
0.427±0.046 
0.426±0.093 


0.552±0.026 
0.557±0.130 
0.562±0.260 


0.415±0.032 
0.569±0.177 
1.040±0.431 


0.081±0.001 0.080±0.001 

0.080±0.001 

0.080±0.001 


0.163±0.018 



this level of precision ijBouchv fc Carrieill2002ll . For the es- 
timation of m CO re we note that the error is introduced at 
two points - first when we determine the scaled core mass, 
(A©/ A) 2 m coro /M, from the calibration curves (Figs. 0] and 
01 and then through the average large separation, A to ex- 
tract m core . The error in A is, however, tiny compared to the 
contribution from the calibration curves while determining 
(A /A) 2 m coro /M. 

For the set of models with M = 1.3Mq,o1 ov = 0.2 we 
have also tested the method for higher error in the frequen- 
cies - relative errors of 5 x 10" 4 and 1 x 10" 3 were added to 
the frequencies. The results, presented in Tabled show that 
the errors are propagated linearly and that the techniques 
to extract m CO ro and X c begin to fail only at the level of 
8v/v ~ 1CP 3 , which is ten times worse than the precision 
level expected from C0R0T. However, in some cases where 
the diagnostic is poorly sensitive to m COTe or X c in a par- 
ticular phase of evolution (e.g., r\ is a poor indicator of X c 
at younger ages) , even a small error in the diagnostic would 
possibly lead to a wrong value of m coro or X c for a star in 
that phase of evolution, even if the errorbar remains small. 
This is what is seen in TableQ for example, in the X c values 
extracted from 77 for relative frequency errors of 5 x 10 -4 . 

3.4-2 Systematic errors 

In applying our technique to determine m coro and X c of a 
given target star, we shall need to use calibration curves 
with approximately the same mass and overshoot as the 



target. In reality, however, these quantities are never a priori 
known exactly. Therefore, we need to investigate the extent 
of systematic errors in m core and X c arising due to the use 
of inexact calibration curves. 

The question of overshoot is one th at asteroseismol- 
ogy hopes to answer. IStraka et al.l i2005f) have suggested 
a method of estimating overshoot through the small sep- 
arations of g-modes. In Section 13.11 we have also demon- 
strated how 77 can be used to detect the presence of core 
overshoot. For this study, however, we realise that the cali- 
bration curves are significantly dependent on overshoot and 
the systematic errors could indeed be large if an inexact 
value of core overshoot is used in the calibration (~ 40% for 
an error of 5d ov ~ 0.2). Therefore, we conclude that an inde- 
pendent estimate of the core overshoot would be necessary 
to use our technique to determine m corc and X c to better 
accuracy. 

For solar-type stars the location on the HR diagram can 
be determined with a fair accuracy which often provides a 
reasonable estimate of the mass of the star. Usually, the 
mass of a solar-type star is known to within O.lOAfQ purely 
from its spectral type and luminosity. We have investigated 
the systematic errors in m corc and X c due to uncertainty 
in the mass by using calibration curves of both 1.2Mq and 
1.3Mq for target models of 1.25AfQ. The results are shown 
in Figs. ll2l and ll3l The errorbars in these figures still indicate 
the random errors due to relative frequency errors of 10 -4 . 
As expected, m CO re is systematically underestimated when 
a lower mass is used for calibration, and vice versa. For an 
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Figure 10. The central hydrogen abundance, as determined 
through the calibration of X c vs n (Fig. [SJ, is plotted against 
the model value. Results are shown for test models both without 
(upper panel) and with (lower panel) overshoot. In each case, the 
calibration curve corresponding to the same mass and overshoot 
was used. The dotted lines show the values which would be ob- 
tained for a perfect match with the model. The errorbars show 
the 3<r error estimates for relative frequency errors of 10 -4 . The 
technique does not work at younger ages for no-overshoot models 
due to the lack of sensitivity of r\ to X c (cf. Fig.|5J. The instances 
of larger errorbars for 1.3 and IAMq occur when the test model 
is close to the kink in the calibration curve. 



Figure 11. The central hydrogen abundance, as determined 
through the calibration of X c vs 9 (Fig. [7J, is plotted against 
the model value. Results are shown for test models both with- 
out (upper panel) and with (lower panel) overshoot. In each case, 
the calibration curve corresponding to the same mass and over- 
shoot was used. The dotted lines show the values which would 
be obtained for a perfect match with the model. The errorbars 
show the 3<r error estimates for relative frequency errors of 10 — 4 . 
The technique does not work at low values of X c due to the non- 
regular variation of 8 with X c at more evolved stages (cf. Fig. 171. 



error of O.O5M0 in mass, the systematic error in m cole /M 
is ~ 0.015 on average. Evidently, the systematic errors in 
m C orc are much larger than the random errors arising due to 
errors in frequencies. 

The situation is, however, different for the determina- 
tion of X c . Fig. 1131 shows that X c can be determined fairly 
accurately from either 9 (only at younger ages) or 77 even if 
the mass is not known exactly. Importantly, there is no sys- 
tematic trend in the error in X c due to an incorrect estimate 
of the mass. The random errors actually span the range of 
uncertainty introduced due to an incorrect estimate of the 
mass. This is not surprising, given the closeness of the cali- 
bration curves for different masses in Figs.© and □ 



4 DISCUSSION 

We have presented the results of our study to relate the 
small frequency separations to the central layers of solar- 
type stars. We have shown that the average small separation, 
(-D02), can be used as a measure of the mass of the convec- 
tive core in stars of mass 1.2-1.4Mq. In addition, we also 
propose two combinations, 9 and rj, of the small separations 
of different pairs of degrees as new diagnostics of the mass 
of the convective core, m corc and the stellar age, in terms of 
X c . We have tested the applicability of our technique to data 
with errors by estimating the random and systematic errors 
in m corc and X c through Monte Carlo simulations. The mo- 
tivation behind the use of such combinations of D02 and D13 
stems from the fact that whil e they both encode slightly di f- 
ferent internal phase shifts jRoxburgh fc Vorontsovll200lfl . 
their combination might be sensitive to the sharp disconti- 
nuity at the boundary of the core. A more thorough the- 
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Figure 12. The mass of the convective core of 1.25Mq mod- 
els, as determined by using the calibration curves of 1.2Mq and 
1.3Mg, is plotted against the model value. The upper panel shows 
results obtained using (-D02), while the lower panel shows results 
obtained using n calibration. The dotted lines indicate the perfect 
match of the model and extracted values. The error bars show the 
random 3<r errors for relative frequency errors of 10 — 4 . The un- 
usually large errorbars correspond to the uncertainty introduced 
by the kink region in the calibration curves (Figs.EJandEJ. 



oretical analysis of these diagnostic quantities is required 
to understand their features better. It turns out that the 
diagnostic value of rj and 6 can be maximised with an ap- 
propriate choice of the frequency ranges indicated here. A 
different choice of the frequency range is not only feasible, 
but might, in fact, be necessary depending on the availabil- 
ity of data. A detailed error analysis using calibration curves 
with the available observed frequencies to be used for the av- 
eraging would, however, be necessary in order to apply the 
technique to real data. 

It turns out that all the calibration curves depend on 
convective overshoot in the core, and therefore we would 
need an independent estimate of the extent of overshoot in 
order to determine m cor c and X c to better than 40%. We find 
that the new diagnostic rj itself is, however, a good indicator 
of the presence of overshoot. The values of rj are much larger 
in the presence of overshoot, irrespective of the mass. 

The mass of the convective core changes as the star 
evolves. For solar-type stars with mass less than 2Mq , it in- 



Figure 13. The central hydrogen abundance of 1.25M0 mod- 
els, as determined by using the calibration curves of I.2M0 and 
1.3M0, is plotted against the model value. The upper panel shows 
results obtained using 8, while the lower panel shows results ob- 
tained using rj calibration. The dotted lines indicate the perfect 
match of the model and extracted values. The errorbars show the 
random 3cr errors for relative frequency errors of 10 — 4 . 



creases with age during the first part of evolution on the 
main sequence before shrinking down slowly in the later 
part. The values of the small separations, however, decrease 
monotonically with age, being more sensitive to the general 
evolution than to the size of the core. We attempt to ac- 
count for the general change in the mean density with age 
by scaling the mass of the convective core with the large 
separation. We find that the average small separation, or 
the new diagnostic rj, varies almost monotonically with the 
scaled core mass, (A0/A) 2 m CO re/M. This scaling enables us 
to calibrate (-D02) or rj against the core mass. It turns out 
that the random error in determining m CO re through (D02) 
is smaller than that through rj. On the other hand, the sys- 
tematic errors due to inaccuracy in the total mass of the 
target star in the calibration curves can be as large as up to 
20% for typical mass errors of 10%. 

We have also shown that the new diagnostic quantities, 
rj and 9 may be used to estimate the age of a main sequence 
star. We find that they serve as nearly complementary indi- 
cators of the central hydrogen abundance - while 6 is a fairly 
good indicator of X c during the early phase of evolution, rj 
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can be used at more evolved stages. Both of these quantities, 
as diagnostics of X c , are fairly insensitive to mass, which 
make them better indicators of the evolutionary stage of a 
star than the small separations themselves. The actual age 
of the star in terms of years would, of course, still depend 
heavily on the stellar mass. The random errors due to errors 
in frequency actually dominate the systematic errors due to 
uncertainty in mass. X c can be determined to 5% for typ- 
ical frequency errors of ~ 0.1/iHz. We note that r\ is more 
sensitive to both m corG and X c when overshoot is included 
in the models. 

While the age of the star determined through r\ or 6 
is quite independent of stellar mass, the determination of 
the mass of the convective core depends on the total mass 
of the star. Further, both of these techniques are sensitive 
to the extent of convective overshoot in the core. We have, 
in this work, provided a method to test the extent of over- 
shoot through the diagnostic r\ itself. The estimate of the 
stellar mass, however, needs to be obtained independently. 
For the fortunate cases of binary stars (e.g.. as in a Cen, 
iPourbaix. Neuforge-Verheecke. fe NoelsUl999T) . the mass is 
usually known to a high accuracy from dynamic consider- 
ations. Even for single stars, a precise location on the HR 
diagram and a good estimate of the metallicity often helps 
to constrain the mass. In any case, a detailed correlation 
analysis between the age, the stellar mass and the extent 
of overshoot will be needed to correctly estimate the uncer- 
tainty in the results obtained from this technique. 
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